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ABSTRACT 

Wc present 12 CO, 13 CO and C 18 J=3^2 maps of the W3 GMC made at the 
James Clerk Maxwell Telescope. We combine these observations with Five Colleges 
Radio Astronomy Observatory CO J=l— >0 data to produce the first map of molecular- 
gas temperatures across a GMC and the most accurate determination of the mass 
distribution in W3 yet obtained. We measure excitation temperatures in the part of 
the cloud dominated by triggered star formation (the High Density Layer, HDL) of 
15-30 K, while in the rest of the cloud, which is relatively unaffected by triggering 
(Low Density Layer, LDL), the excitation temperature is generally less than 12 K. We 
identify a temperature gradient in the HDL which we associate with an age sequence 
in the embedded massive star-forming regions. We measure the mass of the cloud to be 
4.4±0.4xl0 5 M©, in agreement with previous estimates. Existing sub-mm continuum 
data are used to derive the fraction of gas mass in dense clumps as a function of position 
in the cloud. This fraction, which wc interpret as a Clump Formation Efficiency (CFE), 
is significantly enhanced across the HDL, probably due to the triggering. Finally, we 
measure the 3D rms Mach Number, A4, as a function of position and find a correlation 
between M. and the CFE within the HDL only. This correlation is interpreted as due 
to feedback from the newly-formed stars and a change in its slope between the three 
main star-forming regions is construed as another evolutionary effect. We conclude 
that triggering has affected the star- format ion process in the W3 GMC primarily by 
creating additional dense structures that can collapse into stars. Any traces of changes 
in CFE due to additional turbulence have since been overruled by the feedback effects 
of the star-forming process itself. 

Key words: stars: formation - ISM: individual (Westerhout 3) - submillimetre: 
molecular clouds 



1 INTRODUCTION 

One of the most important current questions in the field of 
star formation concerns the effect that environment and, 
especially, feedback may have on the star-formation pro- 
cess, in particular the stellar initial mass function (IMF) 
and the star- formation rate (SFR) or efficiency (SFE). Stars 
appear to form in two main modes. Spontaneous star for- 
mation is the predicted result of the naturally turbulent 
molecular-cloud environment (see e.g. Mac Low & Klessen 
2004; Klessen et al. 2004; Heitsch et al. 2006; Padoan & 
Nordlund 2002) and is expected to produce a low back- 
ground SFR. Triggered star formation, on the other hand, is 
an increase in SFR or SFE due to the effects of a mechanical 
interaction on molecular cloud gas, usually caused by the 
winds, radiation or expanding Hn regions associated with 
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massive stars (see e.g. Elmegreen 1998; Deharveng, Zavagno, 
& Caplan 2005; Ginsburg, Bally, & Williams 2011). 

There are two main ways in which triggering might in- 
crease the SFE locally. The first is by creating new star- 
forming structures (i.e. dense cores), in addition to those 
forming spontaneously in the turbulent molecular cloud gas. 
This mode is most closely described by the collect-arid- 
collapse mechanism (Elmegreen & Lada 1977, Whitworth 
et al. 1994) in which an expanding dense shell, driven into 
a cloud by winds or thermal expansion, becomes gravita- 
tionally unstable and fragments to form dense, star-forming 
clumps. The second mode works by increasing the proba- 
bility that pre-existing dense "cores" will collapse to form 
stars. This would usually require an increase in the ambient 
pressure, either from the passage of a (shock) wave through 
clumpy cloud gas, or when a core is overtaken by an ion- 
isation front (IF). The latter mechanism is described by 
the radiatively-driven implosion (RDI) model (Klein, Sand- 
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Figure 1. MSX image of the W3/4 region at 8 fim. The black con- 
tours trace the 12 CO J=l— >0 outline of the W3 GMC, as observed 
by the FCRAO. The white straight line marks the approximate 
boundary between the HDL and LDL regions. The yellow crosses 
mark the positions of the O stars in the IC 1805 cluster. Several 
regions of interest within the cloud are labeled. 



ford, & Whitaker 1980; Bertoldi 1989; Sugitani et al. 1989; 
Bertoldi & McKee 1990) 

A third potential type of mechanism for influencing the 
SFE would include any that affect the efficiency with which 
dense cores convert mass into stars. This is the same as 
saying that the accretion history of already bound cores is 
affected by a change in the local environment. This might be 
caused by variations in local density and/or effective signal 
speed, which would alter the accretion rate. 

The W3 Giant Molecular Cloud (GMC) offers a poten- 
tial 'laboratory' for constraining the mechanism by which 
feedback triggers new star formation and for quantifying in- 
creases in star-formation efficiency above the spontaneous 
background rate. Both modes of star formation appear to ex- 
ist within W3, each dominating different parts of the cloud. 
Thought to be a prime example of triggered, sequential star 
formation (Lada et al. 1978; Oey & Clarke 2005), it stands 
on the western side of the W4 chimney /sup erbubble whose 
expansion is driven by the winds of the IC 1805 OB asso- 
ciation (Figure 1). This expansion is compressing the east- 
ern side of the W3 GMC and has created a high density 
layer (HDL: Lada et al. 1978) within which there are several 
prominent star-forming regions (i.e. W3 Main, W3 (OH) and 
AFGL333). The rest of the cloud seems so far to have been 
largely unaffected by this interaction and the star formation 
within it should be mainly spontaneous. One notable excep- 
tion is the KR 140 Hn region, located in the far south-west 
corner of the cloud, this may be an example of the spon- 
taneous formation of an isolated massive star which is now 
triggering new star formation in a surrounding shell (Kerton 
et al. 2008). 

Moore et al. (2007) surveyed two thirds of the W3 cloud, 



including all the HDL and the southern half of the remain- 
ing cloud, in the 850-(xm continuum and detected 316 dense 
cores with masses above 13 Mq . Dividing the GMC crudely 
into the two zones, they found that a significantly greater 
fraction of the total gas mass is contained in dense, poten- 
tially star-forming structures in the HDL (25-37%, depend- 
ing on assumptions about the clump mass function, or CMF) 
compared to the diffuse cloud (5-13%), but detected no dif- 
ference in the CMF between the two sections of the cloud. 
These results were interpreted as clear evidence of a collect- 
and-collapse type mechanism at work. However, this result 
was derived assuming a single excitation temperature (30 K) 
everywhere in the molecular gas traced by CO J=l— >0. If 
the gas temperature were significantly higher in the HDL 
than in the remaining cloud, then the contrast in gas mass 
ratio between the two regions may be lower than this anal- 
ysis suggests. 

This paper presents new CO J=3— >2 emission line maps 
of the W3 GMC and an analysis of the physical excitation 
of the cloud molecular gas, in particular the distribution of 
excitation temperatures, using matching CO J=l— >Q data. 
The fraction of gas mass in dense clumps is then estimated 
as a function of position using the existing 850-u.m contin- 
uum results and is compared to the Mach number of the 
turbulence in the CO-traced gas. The paper is structured 
as follows: In Section 2 we detail the data reduction proce- 
dure for the CO J=3— >2 data and describe the CO J=l— >0 
and SCUBA datasets. In Section 3 we describe our analysis 
and discuss the results. Finally, in Section 5, we present the 
conclusions of this study. 

2 OBSERVATIONS AND DATA REDUCTION. 
2.1 HARP Data 

The W3 GMC was mapped with the HARP array receiver 
on the James Clerk Maxwell Telescope (JCMT) on Mauna 
Kea, Hawaii. HARP is a 16-element focal-plane array that 
operates in the 325-375 GHz regime (Buckle et al. 2009). 
The 16 receptors have ~14arcsec beams separated by 30 
arcseconds, giving the instrument a 2-arcminute footprint. 
HARP is combined with the ACSIS digital autocorrelation 
spectrometer backend. 

Observations were made over three consecutive years 
(2006-2008) in 12 CO, 13 CO and C 18 J=3->2 at 345.796, 
330.540 and 329.278 GHz, respectively. All observations 
were taken in good weather with the sky opacity at 225 GHz 
in the range T225 <0.08 (JCMT weather band categories 1 
& 2), using a bandwidth of 250 MHz, giving a basic spectral 
resolution of 26 m s~ . 

As the W3 GMC spans about 1° on the sky, we split 
the cloud into 13 separate tiles of 20 x 20 arcminutes, each 
requiring about one hour to map. All tiles were observed 
using continuous raster scanning and pointing observations 
were made between each tile. We used a sample spacing of 
7.5 arcseconds and the raster scan spacing was half an array 
footprint. Scans were aligned at a position angle of 70° to 
the Declination axis to match the known geometry of the 
cloud. We also observed small raster scan maps of CRL 618, 
CRL 2688 and W3 (OH) to calibrate the science maps. The 
calibration factors applied in the segments of the maps were 
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Figure 2. T^j CO J=3— >2 emission from the W3 GMC, integrated in the range —65 < Vlsr < — 25kms , showing the total area 
surveyed in this transition. Several star-forming regions have been identified on the map. The grid denotes Galactic coordinates. 



between 1.12 - 1.5 for 12 CO, 1.7 - 3.04 for 13 CO and 1.02 - 1.3 
for C ls O. The System Temperature varied between 233 K 
and 283 K for 12 CO with a median value of 242 K; 289 K to 
902 K for 13 CO with a median of 359 K; and 298 K to 340 K 
for C ls O, with a median value of 324 K. The mean pointing 
error was 2.43T0.33 arcseconds for all the observations. 

The observing procedure differed slightly over the years 
as new observing modes became available. In particular, 
the great majority of the 12 CO and 13 CO maps have been 
scanned only along one direction, as the "basket weave" 
mode of orthogonal scanning was not available at the time. 
The C ls O map and later parts of the 13 CO data were 
made using this mode. The velocity range of the cubes is 
-120 kms" 1 to +30kms _1 . 

The raw data cubes were filtered for spikes and, in 



CO, were binned by a factor of nine in the spectral axis 
to achieve a rms noise level of ~0.7K in a 0.23-kms~ wide 
channel. The 13 CO and C ls O J=3->2 maps were binned by 
a factor of 15 to obtain rms noise of ~0.4K in a 0.39-kms _1 
wide channel. The maps were spatially regridded with pixel 
resolution of 7.7 arcseconds and the spectral baselines were 
removed. 



2.2 FCRAO Data 

The 12 CO J=l->0 observations at 115.271 GHz were made 
on 1999 May 7-14 and 2000 April 24-28 at the Five Colleges 
Radio Observatory (FCRAO) 14-m telescope using the 16- 
element SEQUOIA array receiver and FAAS backend with 
256 channels and 80-MHz bandwidth giving 0.8-kms -1 res- 
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Figure 3. The different isotope integrated T R (Kkm/s) emission line maps of the W3 GMC. Top Left: 12 CO J=3->2; Top Right: 
12 CO J=1^0; Middle Left: 13 CO J=3^2; Middle Right: 13 CO J=l->0; Bottom Left: C ls O J=3->2; Bottom Right: C ls O 
J=l->0. 
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olution. System temperatures were in the range 600 K to 
800 K producing rms noise of T\ ~ 0.1K. the rms pointing 
correction was 5 arcseconds or less. These data and their 
reduction are more fully described in Bretherton (2003). 

The 13 CO and C 18 J=1^0 lines at 110.201 and 
109.722 GHz, respectively, were observed simultaneously 
using the expanded 32-element SEQUOIA array on the 
FCRAO 14-m on 2004 March 16-23 March in "On-The- 
Fly" continuous raster mode. The target region was covered 
with 32 individual submaps which are fully Nyquist sam- 
pled on to a 25" grid with a spatial resolution of 49". The 
spectrometer was used with 50-MHz bandwidth centred on 
Vlsr = — 40kms~ 1 , which results in a velocity resolution of 
133 m s _1 . System temperatures were in the range 50 — 80 K 
for the duration of the observations. A fuller description of 
these data can be found in Allsopp (2011). 

All antenna temperatures have been changed to the T£ 
scale using Tp = TJ/^f ss (Kutner & Ulich 1981), where ?7f ss 
is the forward scattering and spillover efficiency, taken as 
0.77 for HARP on JCMT (Buckle et al. 2009) and 0.70 for 
SEQUIOA on FCRAO (Heyer et al. 1998). 



3 RESULTS AND ANALYSIS 
3.1 The W3 GMC 

Emission from the W3 cloud was found in the range —60 < 
Vlsr < —30 kms -1 . The velocity structure and dynamical 
state of the cloud will be discussed in detail elsewhere. Fig- 
ure 2 shows the 12 CO J=3— >2 emission integrated between - 
65 km s _1 and -25 km s _1 . These data cover the whole GMC, 
while the 13 CO and C ls O data in this transition cover 
a more limited area, slightly vignetting the northern and 
southern edges of the cloud. The integrated emission in these 
lines and in the three J=l— > transitions are shown in Fig- 
ure 3. 

Figure 2 shows that, in addition to the warm dense 
gas around the active star-forming regions, 12 CO J=3— >2 
also traces most of the diffuse, extended emission seen in 
lower-level transitions (Figure 3, top right panel). This is 
somewhat surprising, since the critical density of the J=3— >2 
transition is expected to be between 5 x 10 4 cm" 3 at T — 
40 K and 4 x 10 5 cm" 3 at 10 K (Flower & Launay 1985), 
slightly higher than that of CS J=l— >0, and the J=3 energy 
level is E/k = 32.8 K above ground. The transition should 
thus trace the relatively warm, high-density gas associated 
with recent star formation. This may be explained by photon 
trapping caused by high optical depths which may be reduce 
the effective critical density. 

Following Lada et al. (1978) we identify as the high 
density layer the eastern edge of the cloud, adjacent to the 
bubble blown by the IC 1805 OB cluster. The W3 North, 
W3 Main, IC 1795, W3 (OH) and AFGL 333 regions are lo- 
cated here. For the purposes of our analysis, and consistency 
with Moore et al (2007), this is separated from the rest of the 
cloud (which we call the low density layer, LDL and includes 
the HB3, KR140 and the trilobite regions (see figure ??)) 
by a line defined as b = 1.2089 x I + 162.7235, in Galactic 
coordinates. This division is somewhat arbitrary and based 
on the visible extent of the intense star formation in the 
eastern portion of W3. However, while the definition of trig- 
gered and non-triggered cloud regions is not so clear cut, it 



is likely that the feedback effects from the W4 Hn region 
will decrease in strength with distance from the IC 1805 OB 
cluster. Given that the average integrated intensity in the 
HDL is three times higher than that of the LDL (23 K km/s 
compared to 7 K km/s), we assume, in this paper, that the 
expansion of the W4 Hii region has not yet affected these 
regions of W3. Therefore, star formation within the HDL, 
as defined by the equation above, is assumed to be trig- 
gered, whereas the LDL region is assumed to be dominated 
by spontaneous star formation. 

The three brightest and best-known star forming re- 
gions in the cloud (W3 Main, W3 (OH) and AFGL 333) are 
easily identified in Figure 2, running from north to south 
along the eastern edge of the cloud. W3 Main (I = 133.7095°, 
b = 1.2500°) is the most prominent of these and the bright- 
est source in the region at many wavelengths. OH and H2O 
masers have been detected towards sources IRS 4 and 5 
(Gaume & Mutel 1987; Claussen et al. 1984) and it is the 
richest Hn region cluster known within 2 kpc from the Sun 
(Megeath et al. 2008). We find a mean integrated intensity of 
J r TRdV=307.5Kkms~ 1 with a peak of 700Kkms~ 1 in the 
central area and ~55 K kms" 1 in the more diffuse surround- 
ing emission. The cloud associated with W3 Main is seen to 
have rather complex physical structure, possibly shaped by 
the cavity situated immediately to the south which has been 
created by the winds from the young IC 1795 stellar cluster 
(Megeath et al. 2008). 

On the other side of this cavity lies W3 (OH) (I — 
133.9515°, b = 1.0600°), another active high-mass star- 
forming region and also the host of OH and H2O masers that 
point towards two centres of high mass star formation sep- 
arated by 0.07 pc (Norris & Booth 1981; Forster, Welch, & 
Wright 1977). The molecular cloud associated with W3 (OH) 
is clearly elongated north-south, which is consistent with it 
being part of the large-scale, compressed shell produced by 
the winds of the nearby IC 1804 OB association responsible 
for the W4 Hn shell. Its brightest regions have a mean in- 
tegrated 12 CO J=3^ 2 intensity of / T^dv = 167Kkms _1 
with a peak of 300 K km s _1 , while the surrounding structure 
has a mean of ~80Kkms _1 . 

South of W3 (OH) is located the third active star- 
forming region in the HDL, AFGL 333 (Z = 134.2030°, 
b = 0.7630°). This cloud has a less well defined central peak 
than either of W3 Main or W3 (OH) and has a mean inte- 
grated intensity of 49Kkms~ 1 . 

The rest of the W3 GMC contains less intense emission 
and less active star formation. On the south-eastern corner 
of the GMC, the cloud associated with the KR 140 Hn bub- 
ble (Kerton et al. 2008) is easily identifiable. 

Between KR 140 and AFGL 333 is a region we term 
Loops. The CO emission in this area appears diffuse and gen- 
erally has low integrated intensity (J T^dv = 15Kkms _1 ). 
In the 850-(i.m continuum, it appears as a long, fine, looped 
filament (Moore et al. 2007) and in Spitzer data is revealed 
to contain a string of infrared sources (Polychroni et al. in 
preparation). Above KR140 there is a region we call Trilo- 
bite. It has a mean integrated intensity of 14.4 K km s -1 . 
Here, again, Moore et al. (2007) find a number of dense 
cores at 850 urn, indicative of star formation. While this re- 
gion seems rather cut off from the rest of the cloud there is 
a bridge of diffuse material that connects it to the KR 140 
bubble. 
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Figure 4. The 12 CO, 13 CO and C ls O J=3->2 emission lines extracted from the cubes at I = 133.7095°, b = 1.2149° (W3 Main), 
I = 133.9514°, b = 1.0599° (W3 (OH)), I = 133.2028°, b = 1.7628° (AFGL333) and I = 133.4300°, b = 1.4233° (KR140) 



In Figure 4 we plot the emission lines of the three iso- 
topes in J=3— > 2 from individual pixels in the four main 
regions of the cloud (W3 Main, W3(OH), AFGL333 and 
KR140). It is clear from the spectra that the 12 CO J=3^2 
emission line is optically thick and self-absorbed practically 
throughout the cloud. 13 CO J=3— >2 also tends to be opti- 
cally thick and self-absorbed, but only in the brightest re- 
gions like W3 Main or W3 (OH). On the other hand C 18 
J=3— >2 is very weak and we only observe it towards the 
brightest regions of the cloud (W3 Main, W3 (OH) and 
AFGL333). 

Throughout this study we assume that the CO trac- 
ing the molecular gas of the W3 CMC is in a state of local 
thermodynamic equilibrium (LTE), at least in the rotation 
levels J ^ 3. Where the optical depth is high, as is the 
case for 12 CO emission along most lines of sight, the effec- 
tive critical density will be reduced due to photon trapping. 
This means that, in reality, there may be somewhat different 
critical densities for different isotopologues, and the safety 
of the LTE assumption may depend on the rarity of the CO 
species. We ignore this in the following analysis but discuss 
it further below. 

The 3D cubes of all emission lines were collapsed along 
the velocity axis between -65kms _1 and -25kms _1 and 
multiplied by the velocity channel width (Table 1) to pro- 
duce integrated-intensity maps for each species and transi- 
tion. The J=3— >2 data were re-gridded to match the J=l— s>0 
maps so that there was a one-to-one pixel correspondence. 

3.2 Optical Depth 

The measured radiation temperature of a source is given, in 
terms of the excitation temperature T x and optical depth r, 



Table 1. The channel widths of 12 CO, 13 CO, C ls O J=l->0 and 
J=3— >2 spectra. 



Molecule 


Transition 


Channel Width (kms" 1 ) 


12 CO 


l->0 


0.8126 


13 CO 


1^0 


0.1328 


c 18 o 


1^0 


0.1333 


12 co 


3^2 


0.2381 


13 co 


3^2 


0.8301 


c 18 o 


3^2 


0.3348 



by the solution to the equation of radiation transfer, in the 
absence of a background source: 

T R = J(T X ) (l-e- T ) (1) 

where, in LTE, 

J(T:c) = T (eWfcT* _ i _ e WfcT bg n) ' (2) 

v is the frequency and Tb g is the temperature of the cosmic 
microwave background (2.73K). Hence, if the same T% is 
assumed, the ratio of the line brightness temperatures in the 
same transition from two different isotopic species is given 
by 

where r is the optical depth of the more abundant species 
and X is the abundance ratio. 

We adopt a value of X = 77 for the 12 CO/ 13 CO abun- 
dance ratio (Schoier et al. 2002). In this case, the results 
are not very sensitive to the choice of the value of X, par- 
ticularly where the 12 CO optical depth is high. Assuming 
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Figure 5. The optical depth distributions for 12 CO J=3— >-2 (red) 
and J=l— >0 (blue). For the 12 CO J=3— >2 line we measure the 
mean=17, median=14 and the mode=11.2. For 12 CO J=l— >0 wc 
measure the mean=25, median=22 and the mode=10.8. 



Figure 6. The solutions of the excitation temperature, T^, as 
given from the 12 CO J=3— >2 and J=l— >0 opacity ratios (solid 
line) and brightness temperature ratios (dashed line). 



t( 12 CO) 3> 1, the numerator of equation 3 becomes ap- 
proximately equal to unity, providing a first estimate for 
an iterative solution. We use this first estimate and the 
Newton-Raphson iterative method to solve equation 3 to 
calculate the velocity-averaged optical depth per pixel in 
12 CO and 13 CO J=3->2 and J=l->0. We find that the 
range of r( 12 CO) in both transitions is between 5 and 
90 (see Figure 5). At high temperatures we expect that 
-r 3 2( 12 CO)/rio( 12 CO) = 9. However, globally, we find that 
this ratio is much lower, consistent with low gas tempera- 
tures. 13 CO emission is found to be optically thin across the 
cloud in both transitions, with the exception of a few pixels 
in J=3— >2 in the brightest regions (e.g. the central parts of 
W3 Main). 



3.3 Excitation Temperature 

The excitation temperature T x parameterises the relative 
energy-level populations according to the Boltzmann distri- 
bution. Under the LTE assumption, T x is equal to the ther- 
modynamic temperature of the gas. Where optical depths 
are determined for two transitions of the same species, the 
ratio of the two can be used to derive T x , i.e. for the 12 CO 
J=3— >2 and J=l— >0 transitions, 



and J=l— >0 line strengths is given by 



r 32 ( 12 CO) 
rio( 12 CO) 



3 e 



-16.60/Ta- 



-16.60/Ta, 



-5.53/Ts. 



(4) 



This recipe is obtained directly from equation A9, assum- 
ing Tji is measured over the same velocity interval and dis- 
tributed similarly over that interval for both transitions 
j — > i. Equation 4 cannot be solved analytically. Instead, 
we used it to compile a look-up table of optical depth ratio 
values as a function of T x in the range 3 to 34 K (Figure 6). 
The look up table has a resolution of 0.5 K. 

For those pixels in which complete optical depth in- 
formation was not available, usually due to the inadequate 
detection of 13 CO emission, we estimated T x from the ratio 
of line brightness temperatures, using a low optical depth 
approximation, as follows. 

Assuming r« 1, the ratio of the observed 12 CO J=3— >2 



r fl ,32 = 3 r 32 ( 12 co) 

Tr,io 



,16.60/Ta, 



- 1 



- 2.29 x 10" 



tio( 12 CO) ( e 5.53/T x _ iy 



0.152 



(5) 

in which T32( 12 CO)/tio( 12 CO) is given by equation 4. A 
look-up table was again used to obtain T x estimates from the 
observed line brightness ratios. Figure 7 shows the distribu- 
tion of derived excitation temperatures. The spatial resolu- 
tion of the map is that of the J=l— >0 CO data, namely 44". 
The temperature resolution is 1 K, which is the resolution of 
the look-up table, and is less than the uncertainties in the 
data. Where T x > 10 K this error is less than 10%, while in 
regions of lower T x ), the error is around 20%. 

To estimate the error on the excitation temperature we 
calculated the average spectral noise per pixel for each tran- 
sition. These noise maps were added to the integrated in- 
tensity maps and the calculation of r and T x repeated. The 
error estimates quoted above are the average difference be- 
tween the results in the nominal maps and those with added 
noise. 

The majority of the CO-traced gas in the LDL has ex- 
citation temperatures around 8 - 10 K with small excur- 
sions to slightly higher values near star-forming regions (e.g. 
the trilobite and the KR 140 bubble). The three main star- 
forming regions in the HDL all show significantly enhanced 
values of T x in the region of 15 - 30 K. One feature to note 
is that, while the T x distributions in W3 Main and W3 (OH) 
tend to peak centrally, those in AFGL 333 peak near the 
eastern edge of the cloud facing the W4 Hn region (Figure 
7). 



3.4 Gas Mass Distribution 

Having determined T x and r, the column density N is ob- 
tained from equation A8. For the 12 CO J=3— >2 transition, 
and v in kms -1 , the relationship is 



Acq _ 7.67 x 10 17 (T x + 0.922) 
n^ 2 ~ (i^e-ie.eo/T,) 



,16.60/T,, 



t(v)cLv. (6) 
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Figure 7. The distribution of CO excitation temperature in the 
W3 GMC, extimated from the optical depth and integrated in- 
tensity ratios of 12 CO J=3— »2 and J=l— »0. The axes are Galactic 
coordinates. 
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Figure 8. The distribution of molecular gas mass (in Mq units) 
per pixel (0.43 pc/pixel) across the W3 GMC. 

This is converted to molecular hydrogen column density us- 
ing an abundance ratio of [ 12 CO]/[H 2 ] = 9.8 x 1CT 5 (Frerk- 
ing, Langer, & Wilson 1982) and to mass per pixel (Figure 
8) assuming a distance to the cloud of 2kpc (Xu et al. 2006; 
Hachisuka et al. 2006). Integrating over the map, we find 
that the W3 GMC has a total mass of 4.4 ± 0.4 x 10 5 M©, 
consistent with previous estimates of the cloud's mass. The 
error on the mass was estimated in the same way as for T x , 
as described above in Section 3.3. 

The accuracy of the absolute values in Figure 8 is lim- 
ited by the uncertainty in the value of [ 12 CO]/[H2]. The 
value we use is derived from extinction measurements of 
the nearby Taurus and p Ophiuchi molecular cloud (Frerk- 
ing, Langer, & Wilson 1982). Direct measurements of this 
abundance (e.g. Lacy et al. 1994; Watson et al. 1985) can 
differ by factors of 3 to 5 which makes other sources of er- 
ror, propagated from the r and T x calculations, insignificant. 



However, this is a systematic error, and does not affect com- 
parisons made within the map. In terms of random error, 
the most well-behaved regions are those with the highest 
optical depth (r > 20), where uncertainties are of the order 
±10 M© per pixel. In regions with lower excitation temper- 
atures (T x < 10), for which the uncertainties in the optical 
depth produce higher relative errors, the uncertainty on the 
mass per pixel is, on average, ±20 Mq. This translates to 
an error of as low as 4% in pixels with mass higher than ~ 
250 Mq and as high as 30% in pixels with mass lower than 
~ 60 M . 

The spatial distribution of mass shown in Figure 8 
closely follows that of the 12 CO J=3— >2 emission in Fig- 
ure 2, as well as of the 850-um continuum in the HDL and 
southern part of the cloud (Moore et al. 2007). 



4 DISCUSSION 

The CO J=3— >2 emission line data are similar in mapping 
extent to the existing FCRAO CO J=l— >Q data as well as 
the SCUBA continuum observations of the cloud (Moore et 
al. 2007) . The data contain a very large amount of detailed 
information on the physical state of the W3 GMC. This 
paper concentrates on the distribution of gas temperatures 
inferred from LTE excitation temperatures and the results 
of subsequent calculations of the distribution of mass in the 
cloud. The velocity structure and dynamics of the cloud will 
be the subject of a subsequent paper. 

4.1 Gas temperatures 

4-1.1 T x method 

Our excitation temperature results are obtained by assum- 
ing that LTE conditions apply throughout the GMC, i.e. 
that rotational levels J ^ 3 are thermalised, populated ac- 
cording to the Boltzmann distribution which is dependent 
only on temperature. There are several caveats to this as- 
sumption. The first is that the critical density of the 3— >2 
transition is quite high, as mentioned above, and it is quite 
likely that the mean density in the majority of the gas com- 
prising W3 is less than this. Where this is the case, the 
energy-level populations will be determined by the collision 
rate, and so by both the temperature and density of the gas. 
T x will then underestimate the kinetic temperature T of the 
gas and the observed line radiation temperature Tr will be 
less than predicted by LTE. The effect on predicted column 
densities is complex. Whether or not N is under- or over- 
estimated depends at least partly on the gas temperature. 
At low T, an underestimate of T may cause an overestimate 
of N and at high T the reverse may be true. At tempera- 
tures around 30 K small errors in T will have little effect. A 
second point to note is that the effective critical density of 
a transition may be lower where the optical depth is high 
and photon trapping becomes significant. This is likely to 
apply to the 12 CO transitions across most of the cloud and 
to 13 CO along high column density lines of sight and in the 
line centres. This effect may undermine the assumption of 
equal T x for the same transition of different isotopic species 
used to calculate optical depths. 

These issues notwithstanding, our method of deriving 
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T x from the ratio of the optical depths of two different tran- 
sitions of the same species, or from the line radiation tem- 
perature ratio where optical depths are low, has some ad- 
vantages. Firstly, the value of T32/V10 depends on the popu- 
lations of all four energy levels and so the derived value of T x 
represents the distribution of energies relatively well (being 
equivalent to a fit of the Boltzmann distribution) and should 
give a fairly robust estimate of column density, even if T x 
underestimates the real kinetic temperature of the gas. Sec- 
ondly, equation 1 contains, implicitly, the filling factor (i]s) 
of the emission within the telescope beam, i.e. the fraction 
of the beam area filled by the emitting gas. Although ng 
may not be quite the same for different transitions, it will 
be accounted for, to first order, by using Tr and r ratios. 

In other studies, T x has often been estimated from the 
brightness temperature of a single transition, using equation 
1 (e.g. Buckle et al. 2010). This is the best method in the 
absence of data in other transitions, but the LTE assumption 
then models only the relative populations of the upper and 
lower levels of the one transition used (J=3 and J— 2 in this 
case) and does not account for n g, which may be quite small. 

4.1.2 T x results 

The large-scale distribution of T x revealed in Figure 7 con- 
tains few surprises. In general, we see higher temperatures 
(>20K) near regions of active star formation and cooler gas 
(^10 K) elsewhere. We also see two large-scale temperature 
gradients, one running east-west across the whole cloud and 
the other along a north-south axis through the HDL. Along 
the first axis, temperatures range from ~20-30 K in the HDL 
down to ~4-9K in the central and western regions of the 
GMC. It is well known (e.g. Urban & Evans 2007) that in 
regions of high density the gas is in thermal equilibrium 
with the radiatively heated dust, whereas at lower densities 
molecular cooling dominates and the gas cools down through 
molecular transitions. Therefore, this gradient is likely to be 
the result of lower densities to the west of the HDL as well as 
a gradient in the radiation field intensity. The second gradi- 
ent ranges from ~30K in W3 Main to ~10K in AFGL 333. 
Since gas temperatures may be an indicator of the evolution- 
ary stage of star formation within a cloud, this north-south 
trend may imply an age sequence. Such an age gradient has 
also been suggested by Sakai, Oka, & Yamamoto (2005) who 
observed all three regions in atomic carbon. 

In addition to this, there are some detailed differences 
in the T x distribution within the three bright HDL regions. 
In W3 Main, T x peaks clearly in the middle of the associated 
cloud, coincident with the brightest IR and submm sources. 
The embedded YSOs in W3 Main therefore appear to be 
the dominant heating source in the cloud. The surround- 
ing molecular gas may, in fact, be in the process of being 
dispersed by these centrally formed objects. W3 (OH) has 
a lower average excitation temperature of about 20 K. T x is 
also centrally peaked in this source, although less clearly 
than in W3 Main, and there are also high temperatures 
along its western edge, indicating that external heating may 
be important in this cloud. Finally, AFGL 333 exhibits the 
lowest mean excitation temperature (~10K) of these three 
regions and the T x distribution clearly peaks at the eastern 
edge, which is exposed to the radiation from the IC 1805 
cluster (Figure 1). Assuming the embedded YSOs become 



more dominant heating sources with time, these internal T x 
distributions appear to support the idea of an age gradient 
from north to south along the HDL. 



4.2 Mass Distribution 

Obtaining the distribution of T x has allowed us to derive 
the mass distribution of the cloud with much more accuracy 
than previous studies which assume a single temperature. 
Our new estimate of the total mass of the GMC is (4.4 ± 
0.4) x 10 5 M Q , consistent with that of Moore et al. (2007) 
who obtained (3.8±1.1) x 10 5 M Q from 13 CO J=1^0 data, 
assuming T x — 30 K everywhere. This agreement is despite 
most of the cloud having T x < 30 K which, because we are 
below E/k = 33 K, should produce higher mass estimates. 
We find that the mass is almost equally divided between the 
HDL region and the remainder of the cloud (2.23 x 10 5 Mq 
and 2.19 x 10 5 Mq, respectively), even though the latter 
covers almost twice as much projected area as the HDL. 



4-2.1 Clump Formation Efficiency 

We use the existing SCUBA observations of the cloud 
(Moore et al. 2007) along with the masses derived above 
to calculate the fraction of gas in dense, potentially star- 
forming structures as a function of position in the cloud, i.e. 
the clump formation efficiency (CFE). The CFE is a time- 
integrated quantity and can be written as: 



CFE=— — f ™M(t)dt, 

iWdoud J t=0 



(7) 



where M is the rate of formation of dense-core mass from 
the available gas of the cloud. Therefore a high CFE can 
be the result of either a high average dense-core formation 
rate or of a long integration time. We have calculated the 
CFE for the area of W3 surveyed at 850(J.m by Moore et 
al (2007). Sub-millimetre flux densities are converted to gas 
mass using the standard formula, 



M = 



S V D 2 



K„B„(T d ) 



(8) 



where S v is the integrated flux density at 850 um, D is the 
distance to the cloud, k„ is the mass absorption coefficient 
and B„(T) is the Planck function. Td is the dust tempera- 
ture, for which we assume a constant value of 20 K, represen- 
tative of the dust temperatures of dense cores derived from 
infrared SEDs (e.g. Elia et al. 2010; Paradis et al. 2010). 
We assume a constant value of k„=0.01 cm 2 g _1 (Mitchell 
et al. 2001) and so a constant gas-to-dust ratio throughout 
the cloud. The SCUBA map is then regridded to match the 
CO-traced mass data. The division of these two maps gives 
the distribution of the CFE across the cloud (Figure 9). 

As was found for T x , there are two large-scale CFE 
gradients across the GMC. One runs north-south along 
the HDL where W3 Main with CFE~23% and W3 (OH) 
(~20%) exhibit higher efficiency in gas-to-clump conver- 
sion than AFGL 333 (~3%). In comparison, the spontaneous 
star-forming regions in the western section of W3 have very 
low values of CFE, generally below 1%, creating a second, 
east- west gradient. The most obvious interpretation is that 
the high CFE values present in the eastern regions of the 
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Figure 9. The clump formation efficiency (%) as a function of 
position across the W3 GMC. The CFE is given by the ratio of 
the sub-millimetre mass and the CO-traced molecular gas mass, 
calculated where the 850 /im continuum data was available. 



cloud, are the direct result of the interaction between with 
the adjacent expanding W4 Hn superbubble. This is agree- 
ment with the results of Moore et al. (2007) who found that 
26% of the gas has been converted into dense clumps with 
M 13 Mq in the HDL, compared to 5% in the western 
regions of the cloud. They also found that the Core Mass 
Function (CMF) of the cloud does not change significantly 
between the triggered and spontaneous star-forming regions 
of the cloud. There is clear evidence, therefore, that trigger- 
ing due to the W 4 expansion results in the creation of new 
dense structure (clumps/cores) across the affected regions 
that is able to form new stars. The common CMF of the 
two regions suggests that the triggering does not affect the 
star-formation process any further than that, i.e. triggering 
allows the formation of more more massive stars, but it does 
not alter the shape of the cannonical IMF. 

This is in agreement with the collect-and-collapse model 
(Whitworth et al. 1994), where a trigger, such as the winds 
from massive stars, can create a shock that propagates 
through the surrounding medium, sweeping and compress- 
ing the gas adjacent to the forming bubble, creating new 
dense structure that can become gravitationally unstable 
along its surface on long timescales. 

The higher efficiencies measured in W3 Main, compared 
to the W3 (OH) and AFGL333 regions, can also be inter- 
preted as due to different timescales. Assuming that the trig- 
gering from the W4 bubble affects the three regions in a sim- 
ilar way, i.e. it increases the star formation rate to a similar 
degree, then the difference of a factor of three in the clump 
formation efficiency between W3 Main and AFGL 333 star- 
forming regions may be simply due to W3 Main being older 
than AFGL 333. This is consistent with the results of several 
studies that indicate multiple generations of star formation 
in and around W3 Main (e.g. Feigelson & Townsley 2008; 
Rivera-Ingraham et al. 2011). 



4.3 Mach Number 

The turbulent fragmentation model of star formation neatly 
provides a mechanism for the simultaneous support of 
molecular clouds against gravity on large scales and the for- 
mation of dense, star-forming cores in the collisions between 
turbulent flows on small scales, thus naturally explaining 
the low star-formation efficiency (SFE) usually found (Mc- 
Kee & Ostriker 2007). A relation should therefore exist be- 
tween turbulence and the CFE in a given cloud. Padoan 
& Nordlund (2009) use the McKee & Chakrabarti (2005) 
star-formation rate (SFR) model, extended to a magnetised 
medium, to study the relationship between the SFR, the 
virial parameter, a V i r , and the sonic rms Mach number. 
Higher Mach numbers mean stronger shocks, which should 
produce thicker, denser compressed regions as well as addi- 
tional support against gravity on large scales. The relation- 
ship, therefore, is not a simple one. However, the models 
predict a weak negative correlation between the star forma- 
tion rate and the Mach number for turbulence-dominated 
star-forming regions. This implies that there should be more 
large-scale support against gravity where the turbulence is 
stronger, as expected. 

In order to investigate this prediction, we have mea- 
sured the velocity width, a co , of the J=3— >2 emission. Un- 
fortunately, the emission line widths are very much depen- 
dent on whether the line is self-absorbed and/or optically 
thick. Thus, in this analysis we cannot use the 12 CO transi- 
tion since it is largely both self-absorbed and optically thick. 
13 CO also suffers from optical depth effects towards the 
densest and brightest regions of the cloud and we, there- 
fore, have to be very careful when using it to derive the 
velocity widths. Generally, C ls O is the best choice as it is 
optically thin, however, it is a weak line and we only detect 
it in the HDL regions of the cloud. For the 13 CO J=3^2 
we find that the line widths vary between 0.6 and 4.0km/s 
while for the C ls O J=3^2 emission line they vary between 
0.5 and 3.0km/s. 

The total velocity dispersion in the molecular gas can 
be obtained from a co by deconvolving the thermal velocity 
dispersion of the CO molecules, estimated from the sound 
speed ^/3fcT/m co , and convolving the thermal velocity dis- 
persion of the mean molecular gas, i.e. the sound speed 
c a — ^/SkT/^mu- Here [i is the mean molecular weight 
= 2.8, assuming solar neighbourhood abundances. The gas 
temperature T can be estimated using the CO excitation 
temperature T x obtained above (where T x — Tkin assuming 
LTE). This results in c s =0.2- 1.3 kms" 1 through the cloud 
which is systematically lower than the line widths calcu- 
lated above for 13 CO and C ls O, confirming the presence 
of supersonic flows. If the velocity distribution is assumed 
to be Gaussian, we can express the total three-dimensional 
velocity dispersion in the CO-traced gas as 



°3D (total) — 3 



feT 



(9) 



where /i co = 29 and 30 for 13 CO and C 18 0, respectively. 
The Mach number, M, is given by: 



M 2 = cr 3 D( total ) 



(10) 
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Figure 10. The Ml normalised distributions for the HDL region 
as derived from 13 CO (red line) and C ls O J=3— ¥2 (green dashed- 
dot line) and for the LDL region as derived from 13 CO J=3— >2 
(blue dashed line). 



Using equation 9, this can be written 

2 

M 2 = 3.4 x 1CT 4 - x (11) 

where x = 0.903 for 13 CO and 0.907 for C ls O. 

Figure 10 show the range of Mi values in the LDL as 
derived from 13 CO (blue dashed line) and in the HDL from 
13 CO (red line) and C ls O (green dash-dot line). From the 
13 CO it is clear that while there are higher values of M in 
the HDL region the distribution does not differ significantly 
from that of the LDL. This is due to the effects of optical 
depth as well as higher noise in the more diffuse gas that the 
13 CO traces in bot the HDL and LDL regions. The C ls O- 
derived MAch number distribution, on the other hand, does 
not suffer from optical depth effects and suffers less from 
noise levels as it comes from the brightest regions. Therefore 
it should be more indicative of the real Ml distribution in 
the densest regions. 

The first panel of figure 11 shows M plotted against 
CFE for the LDL region from 13 CO J=3->2 (blue stars) 
and the HDL region from both 13 CO (red crosses) and C ls O 
J=3— >2 (green triangles). Note that we have included the Mi 
calculated from all the pixels in all three maps, where there 
is signal, i.e. there are two velocity dispersion estimates for 
regions that have signal in both the 13 CO and C ls O maps. 
Comparisons between the HDL and LDL regions are only 
possible in the 13 CO as we have practically no detections 
in the LDL in C ls O. To determine whether there is a cor- 
relation between Ml and CFE we used the Spearman Rank 
correlation test for all the subsets. For the 13 CO data in 
the HDL we find that despite the spread there is positive 
correlation significant at the 3 sigma level (p=0.25, t=6.7, 
N=675). In the LDL region, on the other hand we find that 
there is no correlation between Ml and the CFE (p=0.066, 
t=1.2, N=240). The spread in the 13 CO J=3->2 derived Mi 
is because of the more diffuse material that surround the 
densest regions in W3 Main, W3 (OH) and AFGL 333 trig- 
gered regions. To minimise this effect and also to negate the 
effect introduced due to the 13 CO lines being optically thick 
in the densest regions, we use the C ls O line, as it only traces 
the most dense regions and is also optically thin. We find a 
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tighter correlation between Ml and the CFE significant at a 
3 sigma level (/9=0.47, t=850, N=261). 

We therefore find no evidence in the current data to 
support the prediction of Padoan & Nordlund (2009). It is 
possible that the lack of correlation between Mi and CFE in 
the western regions of the cloud, which should be dominated 
by spontaneous star formation, could be due to the small 
range in CFE values found there. The significant positive 
correlation found in the HDL can have two interpretations. 
It could be the result of the expansion shock produced by the 
W4 shell interacting with the gas in these regions, increas- 
ing the turbulent velocities and the efficiency with which 
the cloud has formed dense structure, in a way that is in- 
consistent with the predictions of turbulent (spontaneous) 
star formation. A more likely interpretation is that the in- 
creased degree of turbulence is an effect of the higher CFE 
present in these regions, coupled with their age, i.e. it is due 
to feedback from the on-going star formation injecting mo- 
mentum or kinetic energy, which results in higher measured 
Ml. 

In the bottom panel of figure 11 we again plot Ml, as 
derived from the C ls O J=3— >2 emission line against CFE 
per pixel, but separate the data into the three main HDL 
star-forming regions, W3 Main (red crosses), W3 (OH) (blue 
stars) and AFGL 333 (green triangles), revealing that the 
relationship between these quantities is different for each 
region. Spearman rank correlation tests on the three sub- 
samples produce the results p=0.47, 0.51 and 0.65 with t= 
5.6, 4.06, 7.4, and N=115, 63, 77 data points respectively, 
which are all significant at a level > 3a. 

A least-square fit to each of the three regions gives 
different results for each of them with gradient and in- 
tercept values of m=0.11±0.02 and /3=0.55±0.01 for W3 
Main, m=0.17±0.03 and /3=0.56±0.01 for W3 (OH) and 
m=0.26±0.03 and /3=0.57±0.01 for AFGL 333. We note, 
further, that while there is only a small difference of 1- 
2 a between the slopes of W3 Main/W3 (OH) and of 
W3 (OH) /AFGL 333, there is a larger difference of 3-4 a 
between the slopes of W3 Main and AFGL 333. 

This change in the steepness of the slopes in the above 
correlations for the three different regions can be related to 
the evolutionary stage and state of each. If the observed tur- 
bulence were simply due to the mechanical interaction with 
the W4 expansion, and the CFE were determined by the 
turbulence, then all three regions should have more or less 
the same correlation between Ml and the CFE. The change 
in the correlation between the three regions can be inter- 
preted as a result of the on-going star formation affecting 
the cloud. Feedback processes, like outflows or stellar winds 
from forming stars, have been acting on these regions for 
different lengths of time. In W3 Main the slope of the corre- 
lation is flatter than in W3 (OH), which again is less steep 
than that found in AFGL 333. Since W3 Main is thought 
to be more evolved than W3 (OH) and AFGL 333 it follows 
that the flattening of the slope is likely to be because of CFE 
values produced by a longer star-formation timescale. 



5 CONCLUSIONS 

We have presented new CO J=3— >2 maps of the W3 giant 
molecular cloud, obtained using HARP on the JCMT. In 
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Figure 11. Top: M per pixel plotted against the Clump For- 
mation Efficiency in the HDL (red crosses for the 13 CO derived 
M and green triangles for the C ls O derived M) and LDL (blue 
stars for the 13 CO derived M). Bottom: M derived from the 
C 18 J=3— »2 emission line plotted against the CFE across the 
three HDL regions. The fitted lines are the linear least-square 
fits to the three regions. The gradients and intercepts for the 
three regions are: m=0.11±0.02 and /3=0.55±0.01 for W3 Main, 
m=0.17±0.03 and /3=0.56±0.01 for W3 (OH) and m=0.26±0.03 
and /3=0.57±0.01 for AFGL333. The dashed lines indicate the 
1-(T confidence interval in the fit. 



conjunction with FCRAO CO J=l— >Q data, we have used 
these maps to derive the gas properties as a function of 
position within the W3 CMC. 

We have used the ratio of the optical depths of the 
two transitions of 12 CO (where r > 1) and the ratio of 
the brightness temperatures T„ (where r < 1) to derive 
the distribution of excitation temperature in the CO-traced 
molecular gas (Figure 7). We find high excitation tempera- 
tures (T x > 12 K) in the eastern HDL region, where there 
is active star formation. In the remainder of the CMC the 
temperature rarely rises above 10 K. We see a temperature 
gradient along the HDL, where star formation has been trig- 
gered by compression due to expansion of the nearby W4 Hn 
region. We associate this with an age gradient in which W3 
Main is the most evolved of the main star-forming regions, 
followed by W3 (OH) and AFGL 333. 

Using the excitation temperature map, we have ob- 
tained an accurate determination of the distribution of gas 
mass in the W3 CMC (Figure 8). We find that the cloud 
contains 4.4 ± 0.36 x 10 5 M<7>, half of which is located in 



the HDL region. This value is in agreement with previous 
estimates (3.8 ± 1.1 x 1O 5 M ; Moore et al. (2007)). 

We used existing sub-millimetre continuum observa- 
tions of the cloud (Moore et al. 2007) to measure the so- 
called clump- format ion efficiency (CFE) i.e. the fraction of 
molecular gas in the form of dense, potentially star-forming 
structures (clumps), as a function of position in the cloud. 
We find that, in the regions affected by the expanding W4 
Hn superbubble, the CFE has values of 3-25%, much higher 
than that of the rest of the cloud that remains apparently 
unaffected, where values are less than ~1%. We conclude 
that the triggering mechanism that has created the actively 
star-forming HDL in the W3 CMC primarily works by cre- 
ating new dense structures, in agreement with the collect- 
and-collapse model (Whitworth et al. 1994), rather than by 
forcing the collapse of existing structures in the gas. 

We have used the widths of the 13 CO and C ls O J=3^2 
emission lines to derive the sonic rms Mach Number across 
the CMC. We find that there is a positive correlation be- 
tween the Mach Number and the CFE, but only in the gas of 
the HDL. This correlation is opposite to that expected from 
models of turbulence-driven star formation and is probably 
due to feedback from the recent star formation injecting mo- 
mentum into the nearby gas. The slope of this correlation is 
different (Figure 11) in each of the three main star-forming 
region in this part of the cloud (i.e. W3 Main, W3 (OH) and 
AFGL 333) and we interpret this as another indicator of the 
differing evolutionary stages of these three regions. 
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APPENDIX A: DERIVATION OF LTE FORMULAE 

The absorption coefficient of a transition between energy levels j and i at frequency Vji is 

K ( v 3i) = ^ (niBij - UjBji) <j>{vji). (Al) 

Hence, with the relative energy level populations rij/n, populated according to the Boltzmann distribution under LTE and 
with the normal relations between the Einstein coefficients, we have 

= £f>^T (l - e-^ fcT *) 4>{ Vji ). (A2) 



Then, since n,/n is also determined by Boltzmann, 



where Z is the partition function. Now since r(y) = J* k(v) ds ~ k{v)L, where L is the optical path, and = 2j + 1, 



= £(2J + l)^e-^°/^ ^* (l - e - hv »' kT ") H^NiiM (A4) 



Then, integrating over the line, the column density of molecules emitting in the j — > i transition is given by: 



S-7T h "i0 V 2 

N u = — — — — r Z eif# I 1 - e~^H 



J t{v) Av (A5) 



where Z is the partition function and hi>io is the energy of level i above ground. The Einstein A coefficient is given by: 

_ 16tt 3 4 2 j (A) 
where /i is the dipole moment of the molecule. Substituting this in A5 gives: 



The partition function, Z, can be approximated empirically by 



7™ 



Nji = l e ° kC 9 . Z e^V I 1 - e" W ) / T („) di/ (A7) 



Hence, for CO with /i = 0.ff2Debye and replacing the integral over frequency with an integral over velocity v, 

ATji = f^iAi^ ( Ta; + o.922) eT^r 1 - e _ W ) / r(w) du (A9) 



